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ABSTRACT 



PLh We present a study of planet-crossing asteroid (3200) Phaethon at three successive 

PP periheha in 2009, 2010 and 2012, using the NASA STEREO spacecraft. Phaethon is 

,-^ clearly detected in 2009 and 2012, but not in 2010. In both former years, Phaethon 

I brightened unexpectedly by ~1 magnitude at large phase angles, inconsistent with the 

2 ~1 magnitude of steady fading expected from a discrete, macroscopic body over the 

c/j same phase angle range. With a perihelion distance of 0.14 AU and surface tempera- 

I I tures up to ^■lOOO K, a thermal origin of this anomalous brightening is strongly sus- 

pected. However, simple thermal emission from Phaethon is too weak, by a factor >10^, 
J> to explain the brightening. Neither can ice survive on this body, ruling out comet-like 

^^ sublimation. Our preferred explanation is that brightening occurs as a result of dust 

t;^ produced and ejected from Phaethon, perhaps by thermal fracture and/or thermal de- 

^~~J composition of surface minerals when near perihelion. A contribution from prompt 

^T emission by oxygen released by desiccation of surface minerals cannot be excluded. We 

^r^ infer an ejected mass of order 4xl0®amm kg per outburst, where amm is the mean dust 

^^ radius in millimeters. For plausible dust radii, this mass is small compared to the esti- 

K*" mated mass of Phaethon (~2x 10^^ kg) and to the mass of the Geminid stream (10^^ kg 

S^ to 10^^ kg) with which Phaethon is dynamically associated. Perihelion mass-loss events 

like those observed in 2009 and 2012 contribute to, but do not necessarily account for 
the Geminids stream mass. 

Subject headings: asteroids, comets 

1. Introduction 

Object (3200) Phaethon (formerly 1983 TB) is dynamically associated with the Geminid me- 



teor stream, suggesting that it is the long-sought parent of this stream (Whipple 1983). Unlike 



the cometary parents of most streams (e.g. Jenniskens 2008), however, Phaethon has a distinctly 
asteroidal orbit (semimajor axis = 1.271 AU, eccentricity = 0.890, inchnation = 22.2°), with an 
asteroid-hke Tisserand parameter with respect to Jupiter (Kresak 1980) of Tj = 4.5. Phaethon is 
about 5 km in diameter (V-band geometric albedo = 0.17, Veeder et al. 1984) and appears to be 



dynamically associated not just with the Geminids, but also with at least two smaller (kilometer- 
scale) asteroids, namely 2005 UD ( [Ohtsuka et al.||2006t[Jewitt fc Hsiehl|2006t [Kinoshita et al.||2007| ) 
and 1999 YC ( |Kasuga fc Jewitt||2008HOhtsuka et al.||2008D . Together, these objects constitute the 
so-called "Phaethon-Geminid Complex" (PGC) presumably formed by the disruption of a larger 
parent body. Physical observations show that the PGC members share similar neutral-blue colors 
that are relatively rare in the main-belt population and suggesting a common composition. The 
PGC objects all possess small perihelion distances (0.14 AU in the case of Phaethon) resulting in 
high surface temperatures. 

Physical observations of Phaethon when far from perihelion have consistently failed to show ev- 



idence for on-going mass loss, either in gas or in scattered continuum from entrained dust (Cochran 



fc Barker|p84l [Chamberlin et al.||1996t |Hsieh fc Jewittl|2005t [Wiegert et ar]|2008D . Neither do the 
associated bodies 2005 UD and 1999 YC (references above) show evidence for on-going mass loss. 
However, observations of Phaethon at very small solar elongations using the STEREO spacecraft 



revealed anomalous photometric behavior near perihelion in 2009 (Jewitt &: Li 2010, hereafter. 
Paper I). Specifically, Phaethon was observed to brighten with increasing phase angle (near ~80°) 
when near perihelion, a behavior inconsistent with scattering from any macroscopic solid body 



and opposite to the phase functions of known solar system objects (Lane & Irvine 1973 Li et al. 



2007a|b 2009). We interpreted the brightening as caused by an increase in the scattering cross- 



section following the ejection of dust from Phaethon. The required mass of dust is ~2.5xl0 Omm 
kg, where amm is the unmeasured size of the dust grains, expressed in millimeters (Paper I). These 
observations are potentially important both for showing that Phaethon is an active source of matter 
for the Geminid stream, and for illuminating physical processes induced on bodies when close to 
the Sun. 

In this paper, we report new STEREO observations of Phaethon in 2010 and 2012 combined 
with a re-analysis of measurements from 2009 earlier reported in Paper I. The principal question we 
seek to address is whether the anomalous brightening detected in 2009 is recurrent at subsequent 
perihelia. 



2. Observations 



The present observations were made with the Heliospheric Imagers (HI) which are part of 



the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) package (Howard et 



al. 2008; Eyles et al. 2009) onboard the STEREO spacecraft. The HI instruments consist of two 



wide-angle visible-light imagers, HI-1 and HI-2, with field centers offset from the solar center by 
14.0° and 53.7°, and fields of view 20° and 70°, respectively. The HI detectors are charge-coupled 



devices (CCDs) with 2048x2048 pixels. These are usually binned onboard to 1024x1024 pixels, 
resulting in a binned pixel angular size of 70" for HI-1 and 4' for HI-2. The very large pixels subtend 
solid angles 10^ to 10^ times those of pixels commonly used on night-time telescopes, but this is an 
advantage for the intended detection of large scale, diffuse structures in the solar corona. 

We used the standard HI-1 camera Level 1 images in our study. The data are publicly available 
via the UK Solar System Data Center (UKSSDC) web site. Multiple short-exposure images are 
taken before a 1024x1024 pixel-image is transmitted to Earth. The individual exposure time is 
40 seconds, the exposure cadence is 60 seconds and 30 images are combined to make a single 
image having an exposure time of 20 minutes. One such image is obtained every 40 minutes. This 
strategy is designed to remove cosmic rays, achieve statistical accuracy and avoid saturation of the 
background corona. The exposure sequence is chosen so that the drift of the star field through the 
field of view (at ~ 150" hr~^) is comparable to the binned pixel size. The quantum efficiency of 
the CCD camera and the absolute transmission efficiencies of the optics are nearly constant across 
the 6300 to 7300 A wavelength range (Eyles et al. 2009). 



2.1. 2009 and 2012 Observations 

The short orbit period (1.43 yr) offers frequent opportunities to observe Phaethon at perihelion. 
The orbit period of STEREO A is ~346 days so that three STEREO orbits (2.84 yr) are almost 
exactly equal to twice the orbit period of Phaethon (2.86 yr). Therefore, observations in 2009 
(perihelion June 20) and 2012 (perihelion May 02) share similar perspectives viewed from STEREO 
A. 

We obtained the celestial coordinates of Phaethon from NASA's HORIZONS Web-Interface 
and transformed them to pixel coordinates on the HI cameras in the "AZP" Zenithal projection 



(Calabretta & Greisen 2002; Thompson 2006 Thompson, private communication 2012). The sky- 
plane trajectories of Phaethon are shown in Fig. (IT]). As expected, they follow similar paths 
in 2009 and 2012. The angular speed of Phaethon in the images is the result of the combined 
spacecraft parallactic motion and object Keplerian motion. Viewed from STEREO A, Phaethon 
moved relative to the field center at a speed varying between 0" hr~^ and 175" hr~^ from east to 
west, and at a roughly constant speed 180" hr~^ (in 2009) and 200" hr~^ (in 2012) from south to 
north. Phaethon was readily apparent even in a cursory visual examination of the 2009 and 2012 
data. 



2.2. 2010 Observations 

The observing geometry in the intervening orbit in 2010 is completely different from that in 
2009 and 2012. Phaethon was in the field of view of STEREO A from November 11-19 and from 
December 15-25 (upper panel in Figure^. It left the HI-1 field of view six days before its perihelion 



on November 25 and only re-entered the field about twenty days after perihelion in the second time 
period. Phaethon was not detected in either observing window. During the first time interval, the 
phase angle was above 120° making an unfavorable condition for detection. During the second time 
interval, the phase angle was a modest 10°, but the heliocentric distance had increased to ~0.8 AU. 
As a result, the nominal predicted magnitude had fallen to V ~15, making Phaethon too faint to 
be detected in STEREO data. 

Phaethon stayed within the field of the STEREO B HI-1 camera for almost three months 
during its 2010 return (bottom panel in Fig. ^ but still contrived to leave the field of view six days 
before perihelion. Between mid-August and mid-November, the phase angle varied from ~ 0° to a 
maximum 30°, while the heliocentric distance decreased from 1.8 to 0.3 AU. The apparent visual 
magnitude calculated assuming solid body reflection decreased from 18.5 to 12.5. At its predicted 
brightest, Phaethon should have been marginally within reach of STEREO HI-1 camera detection, 
but in practice the object was not detected. 

Presumably as a result of these different observational circumstances, the anomalous bright- 
ening detected in 2009 and 2012 was not found in 2010. Activity specifically at perihelion could 
not be detected because Phaethon was outside the fields of both the A and B cameras when at 
perihelion. We conclude that differences in the observational geometry prohibited the detection of 
Phaethon and its activity at perihelion in 2010. 



2.3. Keck Observations 

To supplement the near-Sun photometry from STEREO, we also observed Phaethon when far 
from the Sun using the Keck 10-m telescope on UT 2012 October 14.5. We used the Low Resolution 
Imaging Spectrometer (LRIS, see Oke et al. 1995) and a broadband R filter (central wavelength Ac 
= 6417A, full-width at haff maximum, FWHM = 1185A) in seeing of FWHM = 1.0". The LRIS 
pixel scale is 0.135"pixel~^. Flat field images were obtained using an illuminated patch on the 
inside of the Keck dome. The data were photometrically calibrated using solar-colored standard 



stars from Landolt (1992). We measured R = 17.25±0.05 at heliocentric and geocentric distances 
2.184 AU and 1.346 AU, respectively, and phase angle a = 18.1°. This measurement refers to an 
unknown rotational phase of Phaethon. 



3. Photometry 

With pixel sizes 70", images of HI-1 are sensitive to the diffuse background coronal emission 
but under-sample the point-spread function of the telescope, and lead to frequent confusion with 
background sources as Phaethon moves across the sky. In these data, the relatively uniformly 
distributed background corona overwhelms both field stars and Phaethon. To suppress the coronal 
background, we grouped images over timespans from a few hours to a few days. Within each 



group, the background corona was relatively constant, and was calculated using the minimum filter 
technique (an IDL code "minjilter" in the SolarSoft IDL package was used) and then subtracted 
from the images in the group. The coronal filtering leaves rapidly varying background structures 
visible in the images, but they are muted in intensity compared to those present in the raw data. 
Sample corona-subtracted images are shown in Figures Q and Q. 

During the 30 minutes required to accumulate a single HI-1 image, Phaethon moved a maxi- 
mum of 1.8 pixels (126") relative to the CCD. We experimented with different photometry apertures 
in order to examine the effects of trailing and background subtraction. In particular, we checked 
that the measured brightness variations are not related to the angular speed of Phaethon and 
therefore to trailing of the images. Very small apertures are sensitive to the trailing, while very 
large ones achieve poor signal-to-noise ratios owing to uncertainty in the large background signal 
from the corona. By trial and error, we chose a photometry box of 5x5 pixels (350" x 350"). The 
sky background was obtained from the median count in 56 sky pixels defined by a box 9x9 pixels 
(630" X 630") wide, surrounding the 5x5 object extraction box. Results obtained using larger pho- 
tometry apertures were consistent with those eventually used, but with a larger uncertainty due to 
noise in the coronal background. To facilitate comparison with measurements from Paper I and to 
make sure that no systematic effects were introduced during our analysis, we elected to completely 
re-reduce the 2009 data-set, as well as those obtained near perihelion in 2010 and 2012. 

The extracted, raw Phaethon photometry from 2009 and 2012 is shown in Figure ([3]). The 
photometry statistics are given in Table M. To be conservative, we set a photometry threshold at 
three times the mean of the sky, as shown by the horizontal lines in Figure ([3]) . The thresholds lead 
to windows in which the photometry is useful lasting for 5.5 days (from June 17.5 to 23) in 2009 
and 3.8 days (from April 30.5 to May 4.3) in 2012. The corresponding phase angles range from 
29° to 130° in 2009, and from 32° to 105° in 2012. In both years, Phaethon shows an apparent 
brightness surge by a factor of ~3 relative to the pre-surge brightness and lasting for approximately 
two days. Each surge is characterized by a sudden rise and fall, with no "plateau" phase in between. 
Phaethon's brightness surges are too large and too long-lived (the photometry aperture crossing 
time is measured in hours while the brightennings last for days) to be caused by passing stars. 
Contaminating stars of sufficient brightness would, in any case, be apparent in visual examination 
of the images but were not seen. In fact, we discarded images in which bright stars interfered 
with the object. Furthermore, the background is constant with respect to time except for small 
excursions of ±10% due to passing stars (blue circles in Figure [s]). For these reasons, we are 
confident that the brightening in 2009 and 2012 is real and associated with Phaethon, not caused 
by background object contamination or sky subtraction errors. 

To photometrically calibrate the measurements of Figure ([3]), we chose field stars near the 
projected path of Phaethon across the CCD. Field stars were chosen to be within 20 pixels (~ 23.3') 
of Phaethon, with magnitudes 8 < F < 10 and of spectral types FGK. Eleven standard stars in 
2009, and fifteen stars in 2012 were available for the Phaethon flux calibration (see Table ^. 
Reference star photometry was obtained in the same manner as for Phaethon. The 6300 A to 



7300 A passband of the STEREO camera is close to astronomical R-band. However, Phaetlion is 
nearly neutral and the reference stars were selected to be similar in color to the Sun. Therefore, 
our measurements of the brightness ratio are equivalent to V-band measurements save for a small, 
color-correction offset <0.1 mag., which is negligible compared to the uncertainties of measurement. 

Figure Q shows the resulting apparent magnitudes of Phaethon as a function of time in the 
time windows of the Phaethon visibility. The times of perihelion UT 2009 June 20 07:12 and UT 
2012 May 2 07:12 are indicated by arrows and marked with the letter "P" . Shown for comparison are 
the magnitudes predicted by the JPL Horizons program, based on an extrapolation of photometry 
obtained at small phase angles. The approximate apparent visual magnitude is calculated for solid 
body reflection at the given heliocentric and Phaethon-STEREO distances. The first observations 
in each year are brighter than the prediction by about 0.5 magnitude. This difference is physically 
insignificant, since the Horizons brightness prediction is based on an assumed (not measured) phase 
function extrapolated to large angles. However, the brightening with time (and phase angle) shortly 
after perihelion in each year is highly significant. Monolithic, macroscopic bodies fade dramatically 
as a result of phase darkening in this range, opposite to the observed brightening. 



4. Discussion 

To correct for the variations in the heliocentric and Phaethon-STEREO distances, Rau and 
Aau, respectively (both expressed in AU), we use the inverse-square law, written 

m{l,l,a) = mobs-5log{Rau^au)- (1) 

Here rriohs is the apparent magnitude and m(l,l,a) is the resulting magnitude which would be 
observed from Rau = ^au = 1 and phase angle a. Equation ([l]) is plotted against a in Figure 
([5]). Since the individual measurements are very scattered (Figure |4]), the plotted curves show the 
running means of 10% of the data. Also plotted in the Figure is m(l, 1, 18.1°) from our Keck 



observation on UT 2012 October 14, and m(l, 1,37.6°) from |Hsieh fc Jewittj ( |2005[ ). The latter 
two measurements are plotted with error bars of ±0.2 mag. to represent uncertainty resulting from 
their being taken at unknown rotational phase. 

The combined ground-based and space-based data show a trend towards fading m{l,l,a) 
upto a ~ 60°. Specifically, Phaethon fades by ~1 mag. from a = 18.1° to a = 60°, giving a 
linear phase coefficient, ~0.024 mag. degree"^, that is unremarkable when compared with other 
asteroids. At larger a, Phaethon in Figure ([5]) shows sudden distance-corrected brightening starting 
at a = 80° in 2009 and a = 65° in 2012. These phase angles correspond to perihelion in each year, 
while maximum brightness is reached at a = 100° in 2009 and a = 80° in 2012, about 0.5 day 
later. Phaethon fades to invisibility at larger phase angles. As noted in Paper I, this brightness 
variation is unexpected for a solid body viewed in scattered light. For such an object, the brightness 



decreases monotonically as the phase angle increases both because a progressively smaller fraction 
of the surface is illuminated and because the scattering efficiency of the surface decreases as the 
scattering angle grows. 

This is emphasized in Figure Ipn, which compares Phaethon with the scaled phase function of 



the Moon, taken from Lane & Irvine (1973). The plotted lunar curve is an average of the phase 
functions at wavelengths 6264 A and 7297 A in Table V of their paper. We show the difference 
between m(l, 1, a) and the scaled magnitude of the Moon in Figure ([6]). The difference plot shows 
that Phaethon's phase dependence is Moon-like before the onset of the anomalous brightening event 
in each year. While the apparent brightness of Phaethon (Figure El) increased by slightly more than 
a magnitude near a = 80° to 100°, the brightening relative to the phase-darkened nucleus (Figure 
p| is a much larger ~2 magnitudes, corresponding to a factor ~6. The figure also shows that the 
apparent, post-peak fading in Figure Q is consistent with the dimming expected from the phase 
function, not necessarily to loss of the scattering cross-section. 



4.1. Mechanisms 

We first consider and reject several mechanisms that might be implicated in the observed 
anomalous brightening at perihelion. 

The intrinsic brightening is too large (2 mag.) and too long-lived (2 days) to be plausibly at- 
tributed to rotational variation of the projected cross-section of the aspherical nucleus of Phaethon. 
Lightcurve observations show a variation <0.4 magnitudes and a rotational period of only ~3.6 



hrs (Krugly et al. 2002). The possibility that the brightening might be caused by a glint (a spec- 
ular refiection from a mirror-like patch on the surface) can likewise be rejected on the basis of the 
longevity of the event. (Separately, the likelihood that the surface of a rocky asteroid could be 
mirror-like in the optical seems remote). 

As noted in Paper I, the brightening of Phaethon cannot be attributed to comet-like processes 
driven by the sublimation of near-surface water ice. This is because surface temperatures on 
Phaethon are far too high for water ice to survive. The deep interior temperature (identified with the 
blackbody temperature of a body moving with Phaethon's orbitally-averaged heliocentric distance) 
is also too high to permit the survival of buried ice (Paper I). Even if it were present, deeply-buried 
ice would be thermally decoupled from the instantaneous solar insolation (the conduction timescale 
across the radius is >10^ yr), leaving no explanation for why the brightening is correlated with 
perihelion. 

The solar wind kinetic energy fiux onto the surface of Phaethon is Eke = pv^/2 (W m~^), 
where p is the solar wind mass density, and v is the solar wind speed. At the Earth's orbit, the solar 
wind number density is about 10^ m~^. Scaled by the inverse-square law to perihelion at 0.14 AU, 
the density is about A^i=5xl0^ m^'^. The solar wind speed varies with time and radius, but is of 
order v = 500 km/s. Substituting p = pninNi, where p = 1 (for protons), and mn = 1.67 x 10^^"^ 



kg, we obtain Eke ~ 0.05 W m~^. This is tiny compared with the solar radiation flux at perihehon, 
Fq/R'J^jj = 70, 000 W m~^. Consequently, we conclude that the solar wind is a negligible source of 
energy and cannot account for the anomalous brightening by impact fluorescence. 

Could part or all of the measured excess optical brightness be thermal emission resulting from 
Phaethon's high surface temperatures when at perihelion? We obtain a rough lower limit to the 
temperature by considering the case of an isothermal, spherical blackbody in equilibrium with 

— 1/2 

sunlight, namely Tbb = 278Rau ■ At perihelion, Rau = 0.14, we find Tbb = 743 K. A practical 
upper limit is given by the sub-solar temperature on a non-rotating body (or a rotating one whose 
rotation axis points at the Sun), namely Tss = \/2 Tbb, giving Tss = 1050 K for Phaethon at 
perihelion. This temperature range is in good agreement with independent estimates (~800 <T< 



1100 K) from Qhtsuka et al. (2009) 



A detailed calculation of the thermal emission from the surface of Phaethon depends on many 
unknowns and is beyond the scope of the present work. Instead, we set a strong upper limit to the 
possible thermally emitted flux density by assuming that the whole surface of Phaethon (not just 
the sub-solar region) is at Tss, the maximum possible surface equilibrium temperature. Under this 
assumption, the flux density in the STEREO bandpass is calculated from 



f =^ 

/A ^2 



rA2 



KB,{T)dX 



AA 



(2) 



in which B\[T) is the Planck function evaluated at T = Tss, fn = 2.5 km is the effective circular 
radius of Phaethon, A is the Phaethon-to-observer distance, the integration is taken over the filter 
transmission from Ai = 6300 A to A2 = 7300 A, and AA = A2 — Ai. 

For comparison with Equation ([2]), we convert the apparent magnitudes, mobs, into flux den- 
sities, /^, using f° = 3.75 x io-i^+"^obs/2.5) [^^.g g-i ^^^^-2 ^-ij (Drilling and Landolt 2000). The 
results are shown in Figure (JTl), where solutions to Equation ([2| (blue) are compared with the 
observations (red) . Evidently, thermal emission even at the peak sub-solar temperature is orders of 
magnitude too small to account for the anomalous brightening of Phaethon observed in our data. 

To see this a different way, we substituted the measured peak flux densities, /^ ~ 10~^^ 
erg s~^ cm~^ A~^ (c.f. Figure M), into the left-hand side of Equation ([2J) and solved for the 
temperature. We find that values T ~ 1650 K (2009) and 1700 K (2012) would be needed for 
thermal emission to account for the anomalous optical brightening. This is far hotter even than the 
sub-solar temperature on the nucleus at perihelion and, therefore, can be discounted as unphysical. 
We conclude that the perihelion brightening of Phaethon is not due to thermal emission from the 
surface. 

The passband of the STEREO camera includes the 6300A and 6363A forbidden lines of oxygen. 
These are "prompt emissions", formed when oxygen atoms are produced in the excited ^D state 
by photodissociation of a parent molecule, for example water (Festou and Feldman 1981). As 



noted earlier, although water ice cannot survive, water might be bound within hydrated minerals 
in Phaethon and released by desiccation. In the case of water, the photodissociation timescale at 
perihelion is about half an hour, so that any water molecules released from Phaethon would be 
destroyed within a single pixel of the STEREO camera. We estimate the flux density produced by 
prompt emission in oxygen, averaged over the passband of the camera, from 



Aoi] 



aQhc 
47rA2AAA 



(3) 



in which a ~ 10% is the fraction of water dissociations leaving oxygen in the excited ^D state, 
Q is the production rate of water molecules, h is Planck's constant, c is the speed of light, A is 
the Phaethon to STEREO distance, A is the wavelength and AA is the filter FWHM, expressed 
in Angstroms. Setting /| = /^, we find that a water production rate Q ~ 10^*^ s~^ (3x10^ kg 
s~^) would be needed to account for the measured excess flux density of Phaethon. Although the 
required rate of production (which is similar to that of comet IP/Halley at perihelion) seems high, 
we cannot rigorously rule out the possibility that some fraction of the excess perihelion emission 
is caused by prompt emission from oxygen. However, the observation that the fading of Phaethon 
after peak brightness follows the phase function of a solid object (Figure [s]) suggests that gas is not 
the dominant cause of the anomalous brightness. 



4.2. Dust 



The remaining alternative is that Phaethon has ejected dust particles with a combined cross- 
section larger than that of the solid nucleus, as earlier concluded in Paper I. In this scenario, the rise 
in brightness in Figure (fsl) then corresponds to the ejection of dust, while the subsequent decline 
in brightness is naturally explained as fading owing to the ever-growing phase dimming, perhaps 
aided by grain sublimation or disintegration. The natural test of this hypothesis would be to search 
for coma scattered by the ejected dust. Unfortunately, as also noted in Paper I, the limited angular 
resolution and high background surface brightness in the STEREO data make the detection of 
resolved coma impossible. 

A temperature-controlled mechanism for the ejection of dust is strongly suggested; the activity 
is observed at the highest (perihelion) temperatures, and is absent in observations of Phaethon taken 
at substantially larger heliocentric distances and lower temperatures. The reflection spectrum 



of Phaethon has been interpreted in terms of thermally modified hydrated minerals (Licandro 



et al. 2007, Ohtsuka et al. 2009, de Leon et al. 2010) while the depletion of sodium in some 



Geminids provides independent evidence for thermal alteration ( Kasuga et al.|2006 ) . The perihelion 
temperatures on Phaethon exceed those needed to break-down phyllosilicates ( Akai|[l992 ) , and are 



sufficient to induce thermal fracture (Paper I, Jewitt 2012). In this sense, thermal disintegration 
and fracture are plausible sources of the anomalous brightening and Phaethon may be accurately 
labeled a "rock comet" (Paper I). An additional requirement is that dust must be cleared from the 



10 



surface in order for these processes to operate. A regolith of fine particles built up in previous orbits 
will inhibit thermal fracture, since small grains are unable to sustain large temperature differences. 
Likewise, surface materials dehydrated by baking in previous perihelion passages must be cleared 
away from the surface in order for dehydration cracking to remain a persistent dust source. 

In both 2009 and 2012, the apparent brightness increased by ~2 magnitudes relative to the 
nominal phase curve (Figures ^ and ^ , corresponding to a factor of ~6. Given that the cross- 
section of the nucleus of Phaethon is C„ = vrr^ ~20 km^, the cross-section of added dust is then C 
= 100 km^, in both years. The mass in spherical particles of mean radius a having cross-section, 
C, is 

M, ~ ^. (4) 

where p is the grain density. With p = 3000 kg m~^, we find M^ ~ 4x10® cLmm kg) where amm is 
the grain radius expressed in millimeters. The mass of the nucleus, represented as a 2.5 km radius 
sphere of the same density, is a much larger 2x10^^ kg. 



The Geminid stream mass is 10^^ < Ms < 10^^ kg ( Hughes fc McBride|l989 Jenniskens|l994 ), 



while the Geminid stream lifetime estimated on dynamical grounds is r ~1000 year (Gustafson 1989 



Ryabova||2007 ) . Accordingly, to explain the entire mass in the Geminid stream through events like 



those observed here would require a number of similar events per orbit, A^, given by 

where P = 1.4 year is the orbital period. Substituting, we obtain A^ ~ 4a~}„ to 40a~|„. If the 
particles are millimeter-sized, amm = 1, then the Geminids could be supplied in steady-state by 4 
< A^ < 40 outbursts like the one observed, each orbit. 

However, there is no compelling physical reason to assume that Geminid stream production is 
in steady state. While amm = 1 rnay approximately represent the radii of the Geminid meteors, 
much larger examples up to 5 kg in mass (equivalent radius ~7 cm) have been inferred from 



Lunar night-side impact flashes (Yanagisawa et al. 2008). Moreover, Phaethon is but part of the 



Phaethon-Geminid Complex which includes at least two other kilometer-scale bodies caused by 
fragmentation on an all-together much larger and longer (10^ yr?) scale. We conclude that, while 
continuing mass-loss near perihelion may contribute to actively replenishing the smaller Geminids, 
the PGC complex as a whole is likely the product of a more ancient and catastrophic breakup 
(c.f. Jewitt and Hsieh 2006, Kasuga and Jewitt 2008). 



11 



4.3. The Future 

Many puzzles remain in understanding the anomalous brightening of Phaethon and its relation 
to the ejection of dust and to the Geminids. We list the following key questions: 

1. What is the value of the effective dust radius, Cmm? This radius directly affects estimates 
of the ejected mass through Equation Q and so determines the extent to which on- going 
activity in Phaethon contributes to the Geminid stream. 

2. What is the origin of the ~0.5 day phase lag between perihelion and the anomalous brightening 
in both 2009 and 2012? Is the brightening always lagged relative to perihelion by this amount? 

3. How many brightening events occur per orbit and are these events always of the same am- 
plitude? The nature of the STEREO data curtails our ability to detect brightening events 
outside a limited window of accessibility. 

4. Can evidence for mass loss be detected when Phaethon is far from perihelion, despite past, 
failed attempts? Thermal fracture and mineral decomposition are likely only at the extreme 
temperatures found near perihelion. However, slow-moving dust may linger in the post- 
perihelion months, making these the prime time for future attempted dust observations. 

5. Can spectra be obtained at perihelion in order to search for the forbidden lines of oxygen? 

6. What is the mineralogical composition of Phaethon and are hydrated silicates present on its 
surface? 

7. Can meteoroids ejected during recent perihelion events be detected and distinguished from 



older Geminids? Ryabova (2012) reports that this will be difficult. 
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5. Summary 

We report new observations of planet-crossing asteroid and Geniinid meteoroid parent (3200) 
Phaethon, using the NASA STEREO solar spacecraft. We find that 

1. (3200) Phaethon exhibited anomalous brightening when at perihelion in 2009 and 2012, 
but not in 2010 (the latter likely owing to unfavorable observing geometry). The distance- 
corrected apparent brightness increased near phase angle 100° in 2009 and 80° in 2012, in 
both years ~0.5 day after perihelion passage. This brightening lies in stark contrast to the 
monotonic fading expected from phase darkening on a macroscopic body. 

2. The most direct interpretation is that Phaethon brightens because of a sudden increase in 
the scattering cross-section due to the ejection of dust with a mass M^ ^ AxlO^amm kg, 
where Omm is the effective dust radius in millimeters. A contribution from prompt emission 
by atomic oxygen cannot be excluded. 

3. Thermal fracture and the decomposition of hydrated silicates are two plausible mechanisms 
of dust production at the ~1000 K surface temperatures attained near perihelion. Both are 
difficult to quantify in the absence of more detailed information about the composition of 
Phaethon. 

4. Phaethon has only very limited visibility in the STEREO field of view (typically <5 days 
per orbit). The detection of anomalous brightening twice in two favorable observing windows 
suggests that this phenomenon is common. 
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Table 1. Photometry Statistics and Detection Characteristics 



Year Target Mean Std. Deviation Threshold^ 



Datc^^ 



Phase Angle^ 



2009 Sky 0.09 

Phaethon 0.44 



0.04 
0.37 



0.27 



June 17.5-23.0 



29°-130 



2012 Sky 0.10 

Phaethon 0.36 



0.12 
0.40 



0.30 April 30.5 - May 4.3 32°-105° 



TThe adopted threshold for Phaethon photometry, equal to three times the sky mean. 
■fTime interval during which Phaethon's light curves are above the threshold. 
^ Phase angle ranges corresponding to these date ranges. 
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Table 2. Standard Stars 



HlC 



HD SAO [R.A., Dec] Magnitude Spectra 



2009 



1 


47142 


83099 


98670 


144.0 


10.6 


8.50 


F2 


2 


48246 


85143 


117932 


147.5 


9.7 


8.48 


G5 


3 


48381 




98815 


147.9 


9.8 


8.40 


KO 


4 


48624 




98843 


148.7 


9.9 


8.40 


F8 


5 


48866 


86340 


118022 


149.5 


9.7 


8.50 


F5 


6 


49194 




98915 


150.6 


9.9 


8.97 


KO 


7 


50146 


88724 


99018 


153.6 


10.5 


9.40 


G 


8 


50110 


88680 


99012 


153.5 


10.8 


8.02 


05 


9 


50424 


89209 


99041 


154.4 


11.3 


8.60 


F5 


10 


50467 






154.6 


11.8 


9.30 


KO 


11 


50911 


90042 


99096 


155.9 


14.8 


8.58 


05 


2012 














1 


52259 


92458 


118427 


160.2 


5.5 


8.90 


05 


2 


53032 


93981 


118531 


162.8 


5.0 


8.68 


K5 


3 


53165 


94220 


118554 


163.1 


5.0 


8.17 


KO 


4 


53899 


95529 


118622 


165.4 


5.1 


8.80 


05 


5 


54120 




118647 


166.1 


5.1 


9.10 


F8 


6 


54213 




118652 


166.4 


4.9 


8.70 


KO 


7 


54188 




118649 


166.3 


5.2 


9.61 


F8V 


8 


54481 


96762 


118684 


167.2 


5.0 


8.83 


FO 


9 


54660 




118707 


167.8 


5.1 


9.30 


F8 


10 


54788 


97459 


118724 


168.2 


5.2 


8.50 


OO 


11 


55323 


98451 


118791 


169.9 


6.2 


9.20 


KO 


12 


56338 


100362 


118913 


173.2 


9.4 


8.38 


F8 


13 


56542 


100727 


118935 


173.9 


9.6 


8.35 


F2 


14 


56491 






173.7 


10.0 


9.70 


F2 


15 


56646 


100904 


99697 


174.2 


10.5 


8.26 


F5 
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.Ji^l^ 






2012May 



Fig. 1. — The path of Phaethon (white circles) across the field of view of the STEREO A HI-1 
camera in 2009 (left) and 2012 (right). Numbers along the path show the day of the month. The 
perihelia are indicated with the letter "P" and a large circle. The sun is on the right. The images 
were taken on the dates when Phaethon was at the perihelion. Both panels show 400x700 pixels 
(7.8° X 13.6°). 
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Fig. 2. — The path of Phaethon (white circles) across the STEREO A (upper) and B (lower) HI-1 
fields of view in 2010. Numbers along the path show the date in month/day format. Perihelion 
occurred on UT 2010 November 25 18:00, at which time Phaethon was not within the field of either 
STEREO camera. Panels show the full size HI-1 images of 20°x20°. The sun is on the right in the 
upper panel, and on the left in the lower panel. 
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Fig. 3. — Raw photometry of Phaethon (red circles) as a function of time near perihelion in 2009 
and 2012. The median sky brightness surrounding Phaethon is shown with blue circles. Smoothed 
curves (red and blue lines) have been plotted to guide the eye. Horizontal lines represent the 
photometry thresholds that are three times of the sky mean (see Table IT]). Above the levels, 
Phaethon was detected. 
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Fig. 4. — Apparent V magnitude of Phaethon (red circles) in 2009 and 2012. The magnitudes 
predicted by NASA's ephemeris software are shown for reference (blue lines). Letter "P" and 
arrows mark times of perihelion, date = 20.3 for June 2009 and date = 02.3 for May 2012. The 
red curve is a smoothed fit to the data added to guide the eye. The time ranges correspond to the 
valid Phaethon photometry measurements (see Table ([l])). 
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Fig. 5. — Magnitude at i? = A = 1 AU, vs. phase angle (a) in 2009 and 2012. The reduced 
Phaethon magnitudes are plotted in red (2009) and blue (2012) curves, and are smoothed fits to 
the actual data points. Two Keck data points are from a new measurement in 2012 for m(l, 1,18.1°); 
and from Hsieh & Jewitt (2005) for m(l,l,37.6°). The lunar phase function is over-plotted from 



Lane & Irvine ( 1973[ ) with the thick solid black curve. For a > 120°, the lunar phase function is 
extrapolated (dashed curve). The letters "P" indicate phase angles a = 79° (red) and a = 66° 
(blue) corresponding to the perihelia in 2009 and 2012, respectively. Note that these perihelia 
correspond to the starts of the Phaethon brightnesses. 
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Fig. 6. — Phase angle dependence normalized to the phase function of the Moon. The phase angles 
at perihelion are marked with arrows and the letter "P" . 
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Fig. 7. — Flux density as a function of time from photometry (red circles) and the thermal emission 
(blue curves). The latter is calculated from Equation (pi) in 2009 (top) and 2012 (bottom). The 
perihelia are marked by the letter "P" with arrows. 



